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1.  INTRODUCTION 


In  Cll,  Jarvis  at.  al.  praaant  a  aathod  for  daployaant  planning 
eallad  Syataa  lor  Cloaura  Optiaization  and  Planning  < SCOPE).  This 
aathod  usas  Bandars'  dacoaposition  aathod  to  dacoaposa  tha  problaa 
into  aanagabla  ooaponants.  This  raport  givas  tvo  possibla 
dacoaposition  aathods  lor  a  ralatad  ailitary  daployaant  problaa  - 
Intra-CONUS  aovaaant. 

SCOPE  is  principally  concarnad  with  tha  Aovaaant  aovaaant 
raquiraaants  Iroa  ports  ol  aabarkation  (POEs)  to  ports  ol  dabarkation 
(PODs),  using  assats  such  as  airplanas  and  ships.  Givan  a  solution 
to  this  problaa,  it  is  still  nacessary  to  plan  tha  aovaaant  ol  tha 
lorcas  to  thair  raspactiva  POEs.  This  problaa  is  tha  Intra-CONUS 
(lor  CONtinantal  Unitad  Statas)  Traval  Problaa  (ICTP). 

Ona  rasult  ol  SCOPE  is  an  aaaignaant  ol  lorcas  to  POEs.  For 
ICTP,  it  will  ba  assuaad  that  aach  lorca  is  assignad  to  axactly  ona 
POE  (il  SCOPE  splits  a  aovaaant  raqulraaant  aaong  aultipla  POEs  than 
ICTP  vill  hava  aora  than  ona  lorca  assodatad  with  that  aovaaant 
raqulraaant ) .  Each  lorca  is  also  associatad  with  a  uniqua  origin. 

Tha  aovaaant  ol  tha  lorca  vill  ba  Iroa  tha  origin,  through  a  road  or 
rail  natvork,  to  tha  POE.  It  is  assuaad  that  tha  road  and  rail 
natvorks  ara  assantially  uncapncitatad  so  tha  only  constraints  ara  on 


2.  MATHEMATICAL  MODEL 


Daf in*  Xi • .  to  ba  tha  tonnaga  of  f ore*  i  that  arrivaa  at  its  POE 
at  tlaa  t  aaploying  aoda  a,  vhara  a  la  althar  M  (for  aotor)  or  R  (for 
rail).  Let  r(i,  a,  t)  ba  th*  tlaa  that  forca  1  arrivaa  at  Ita  POE  If 
It  laavaa  ita  origin  at  tlaa  t  uaing  aoda  a.  Thia  function  la  alaply 
j  a  ahlft  function  that  aodala  tha  tlaa  naadad  to  aova  froa  tha  origin 
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to  tha  POE. 

EXAMPLE:  If  it  takaa  forca  i  tvo  daya  to  aova  froa  ita  origin  to 
ita  POE  by  rail  than 

r(l,  R,  4)  ■  6 
r(i, R, 5)  -  7 

Lat 

0(1)  ba  tha  origin  for  forca  1) 

E ( i )  ba  tha  POE  for  forca  i ; 

HRi  ba  tha  tonnaga  of  forca  i; 

POEj  •  ba  tha  capability  (in  tona/day)  of  POE  J  at  tlaa  t; 
and  0Rk«a  ba  tha  capability  (in  tona/day)  of  origin  k  at  tiaa  t 
with  aoda  a. 

Tha  ICTP  can  ba  aodallad  vith  tha  following  conatrainta: 

(1)  Ei  £■  Xi«.  2  MR(  for  aach  forca  i; 

(2)  E|  inn./  E.  X| « ,  &  POEj «  for  aach  POE  j  and  tiaa  t; 

(3)  Eniiit.il  Xi,M|l,na  S  ORata  for  aach  origin  k,  tiaa  t  and 
aoda  a. 

(Tha  notation  Emui.j  aaana:  aua  ovar  all  forcaa  i  vith  E(i) 

■ 

Conatraint  (1)  fixaa  tha  aaount  aovad  of  aach  forca  to  ba  tha 
forca  aiza.  Conatraint  (2)  aodala  tha  throughput  capability  of  aach 
POE  during  aach  tlaa  parlod.  Conatraint  (3)  aodala  tha  capacity  of 
aach  origin  during  aach  tiaa  parlod,  vith  aaparata  conatrainta  for 
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rail  and  not or  aovaaant 


Thla  aodal  can  ba  aolvad  aa  a  llnaar  prograa.  Unfortunately,  tha 
aodal  can  ba  quite  large.  A  daployaant  with  10,000  forcaa,  SO 
orlgina,  SO  POEa  and  a  tiaa  horizon  of  100  daya  will  have  2, 000, 000 
variablaa  and  25, 000  conatrainta.  Tha  aodal  ia  far  too  large  to  ba 
aolvad  by  a  general  linear  prograaaing  coda  in  a  raaaonabla  aaount  of 
tiaa.  Tha  following  aectiona  develop  tvo  dacoapoaition  aathoda  which 
exploit  tha  aabaddad  atructura  of  tha  aodal  and  aolva  auch  eaaller 
problaaa. 
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3.  TIME  /  NODE  DECOMPOSITION 


Define  y,  •  to  b*  th*  aRount  of  fore*  i  that  arriv**  at  lta  POE  at 
tla*  t.  Th*n  y, »  ■  x, .  «  ♦  x,**.  Th*  aod*l  becomes! 

( 1 ' >  E,  yt •  2  NR,  for  *aeh  fore*  1} 

(2')  E,  umij  y,»  £  POE,,  for  *aeh  POE  j  and  tin*  t j 

< 3 ' )  E, !•«,,. k  X|Fii,,,«i.  £  ORk»a  for  *ach  origin  k,  tia*  t,  *od*  a 
( 4 ' )  x,,t  ♦  x,  • «  ■  y, ,  for  *ach  fore*  i  and  tia*  t. 

Thia  can  b*  aolv*d  with  B*nd*r*'  d*coapo*ition  a*thod.  Th* 
aa*t*r  probl**  ar*  constraint*  1*  and  2'.  This  probl**  is  a 
transportation  n*tvork  flow  probl**.  Th*  sourc*  nod*s  of  th* 
transportation  probl**  ar*  th*  fore**.  Th*  destination  nod**  ar*  th* 
tia*  expanded  POE*.  A  fore*  is  connected  to  a  tia*  expanded  POE  nod* 
if  it  us*s  that  POE  and  it  can  arrive  at  that  tia*.  Hence,  th* 
network  consists  of  aany  smaller  networks,  on*  for  each  POE  (s** 
Figure  1). 

Given  a  solution  to  th*  aast*r  problem,  y,»',  th*  subproblem  is: 

(  3  )  E  X,  ,  (  |  i  ■  «  t  ,  a  £  ORk  «  ■  } 

(4*  )  X,  ,  m  ♦  X,  «  m  •  y,  a  '. 

This  is  also  a  transportation  network  flow  problem.  Th*  source 
nodes  ar*  th*  tim*  *xpand*cf  movement  requirements  and  th*  destination 

nodes  are  the  time  expanded  origin  nodes,  one  set  for  each  mode.  A 

tia*  expanded  force  is  connected  to  exactly  two  origin  nodes t  on* 
for  th*  corresponding  time  expanded  rail  origin  and  on*  for  th* 
corresponding  tim*  expanded  motor  origin.  In  other  worde,  if  a  force 
node  corresponds  to  arriving  at  its  POE  on  day  5  and  rail  travel  is 
two  days  and  motor  travel  is  on*  day,  then  th*  day  S  force  nod*  will 
be  connected  to  the  day  3  rail  nod*  and  the  day  4  motor  node  (see 
Figure  2).  This  network  also  consists  of  many  smaller  networks,  on* 
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for  •■oh  origin.  Not*  that  many  of  the  fore*  nodes  will  have  aupply 


zero  <for  the  aupply  corresponds  to  y» » # )  further  aiaplifying  th* 
network. 

Th*  solution  to  th*  aubprobl**  g*n*rat**  a  new  constraint  which 
is  will  b*  add*d  to  th*  sast*r  problem.  Iteration  b*tw**n  th* 
probl*as  continues  until  a  "good"  answer  is  generated. 

A  simple  Interpretation  is  available  for  this  decomposition.  The 
master  problem  decides  how  many  tons  of  each  force  will  arrive  at  its 
POE  on  each  day.  Ignoring  the  origin  constraints.  The  subproblem 
uses  this  result  and  decides  whether  to  use  rail  or  motor  movement  to 
•void  congestion  at  the  origin.  This  decomposition  will  be  called 
the  TINE  /  MODE  decomposition. 


4.  CHANNEL/ ASSIGNMENT  DECOMPOSITION 


An  alternative  decomposition  strategy  is  to  create  channels 
between  origins  and  POEa  and  then  to  allocate  forces  to  the  channels. 
This  decosposition  assumes  that  the  transportation  tiae  is  dependent 
on  only  the  node,  POE  and  origin  and  is  independent  of  the  force 
involved. 

Define  yk j , .  to  be  the  tonnage  moved  from  origin  k  to  POE  j  at 
time  t  via  mode  m.  Let  q(k, j,m,  t)  be  the  time  a  force  will  leave 
origin  k  if  it  arrives  at  POE  j  at  time  t  via  mode  m.  This  gives  the 
following  model : 

( 1  ■  >  E»  E.  Xu.  2  MRi  for  each  force  i; 

<2")  Liimi.j  kii.o  x, ».  *  y»  j » ■  for  each  POE  j,  origin  k, 

time  t  and  mode  a; 

< 3* )  Ek  E.  ykj*m  S  POEj «  for  each  POE  j  and  time  t; 

(4")  Ej  yt j , 1 1 , j , ■ , 1 1 *  S  0Rk • >  for  each  origin  k,  time  t  and  mode  m. 

Again,  this  model  can  be  solved  with  Benders'  decomposition 
method.  The  master  problem  consists  of  constraints  3 ■  and  4”. 
Constraint  3”  forces  the  amount  arriving  at  each  POE  during  each  time 
period  to  be  less  than  the  capacity  of  the  POE.  Constraint  4"  limits 
the  amount  leaving  an  origin  to  the  origins  capacity. 

This  problem  is  a  transportation  network  flow  problem;  the 
origins  are  the  sources  and  the  POEa  are  the  sinks.  Unlike  the  TIME 
/  MODE  decomposition  this  problem  does  not  break  into  smaller 
problems.  (see  Figure  3). 

The  subproblem  consists  of  constraints  ia  and  2”.  Given  a 
solution  to  the  master  problem  <ykj«a>,  the  subproblem  attempts  to 
allocate  the  forces  to  the  created  channels  (origin  -  POE  pairs). 

Aga  in,  this  subproblem  is  a  transportation  network  flow  problem;  the 


forcaa  ara  aoureaa  and  tha  ohannala  ara  ainks.  Thia  problaa  conaiata 
of  nuaaroua  aubnatvorka.  Tha  aourca  nodaa  in  aach  aubproblaa  conaiat 
of  thoaa  forcaa  with  a  coaaon  origin  and  POE.  Tha  aink  nodaa 
corraapond  to  channala  opanad  in  tha  aaatar  problaa.  In  ganaral,  tha 
aubnatworka  will  ba  vary  aaall  <aaa  Figura  4). 

A  aolutlon  to  tha  aubproblaa  will  craata  a  conatralnt  in  tha 
■aatar  problaa  and  tha  procaaa  itarataa.  Thia  dacoapoaition  will  ba 
callad  tha  CHANNEL  /  ASSIGNMENT  dacoapoaition  ainca  it  craataa  channala 
and  than  aaaigna  forcaa  to  thaa. 


Origin  1  and  POE  1  channels 
opened  by  Master  Problem 


Origin  2  and  pOE  1 


Origin  1  and  POE  2 


Origin  2  and  pOE  2 


Figure  4.  Suboroblem  for  CHANNEL  /  ASSIGNMENT  Decomposition 
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5.  COMPARISON  BETWEEN  DECOMPOSITIONS  ! 

Both  decompositions  aaparata  ICTP  into  two  or  aora  tranaportation 
network  problems,  which  are  vary  efficiently  aolvabla.  Tha  aizaa  of 
tha  problaaa  ara  vary  dif farant.  Conaidar  a  problaa  with  10, 000 
forces,  SO  POEa,  SO  origina  and  a  tiaa  horizon  of  100  daya. 

Tha  aaatar  problaa  of  tha  TIME  /  MODE  dacoapoaition  will  hava 
IS,  000  nodaa  and  1,000,000  area.  Although  thia  problaa  oonalata  of 
SO  aubnatworka  of  roughly  aqual  aiza,  tha  conatrainta  ganaratad  by  j 

Bandars'  dacoapoaition  aathod  will  link  thasa  networks,  pravanting 
Individual  solution.  Tha  aaatar  problaa  of  tha  CHANNEL  /  ASSIGNMENT 
dacoapoaition  has  IS,  000  nodaa.  Tha  nuabar  of  arcs  is  dapandant  on  I 

further  details  of  tha  problaa  <how  aany  origin-POE  pairs  hava  at 
least  one  force ) ,  but  la  certainly  no  aora  than  500, 000  and  aay  be 
auch  less. 

Tha  subproblaa  of  tha  TIME  /  MODE  dacoapoaition  has  SO  networks 
with  roughly  300  nodes  and  1,000  arcs  each.  Tha  subproblaa  for  tha 
CHANNEL  /  ASSIGNMENT  dacoapoaition  consists  of  aany  vary  saall 
networks.  Tha  aiza  and  nuabar  of  thaaa  networks  la  dapandant  on 
problem  data,  but  will  probably  be  hundreds  or  thousands  of  networks 
with  roughly  tan  nodaa  and  fifty  arcs. 

Both  aaatar  problaaa  ara  vary  large,  and  tha  presence  of  tha 
Bandars'  conatrainta  requires  tha  use  of  a  network  flow  with  aide 
constraints  aathod.  Thia  aathod  ia  slower  than  a  standard  network 
flow  aathod.  Fortunately,  tha  number  of  aide  constraints  will  be 
vary  small  (one  par  iteration)  so  thasa  problems  ara  still  aolvabla. 

Further  assumptions  on  tha  problem  (Ilka  hard  windows  for  tha  forces) 
will  further  reduce  tha  problem  sizes. 


Without  extnesl 


tasting,  it  is  not  possible  to  determine  which 


of  the  two  decomposition*  vlll  provide  the  better  anever  for  a  small 
number  of  Iterations.  Both  methods  seem  to  create  reasonable  sized 
subproblems,  so  both  are  possibilities  for  further  testing. 


6.  USING  THE  MODES  SOLUTION 


This  section  discusses  various  ways  in  which  MODES  solution 
information  can  be  used  as  input  to  the  Intra-Conus  aodels  presented 
in  this  report.  The  POE  assignment,  generated  by  MODES,  is  required 
inforsation  which  MODES  would  have  to  provide  to  either  Intra-Conus 
model.  The  remainder  of  the  information  discussed  in  this  section  is 
desired  information,  generated  by  the  MODES  model,  which  could  be 
used  in  a  mathematical  programming  based  decomposition  model  or  in 
other  model  types  <such  as  a  deterministic  simulation  approach). 

The  solution  generated  by  MODES  allocates  forces  to  specific  time 
expanded  channels.  The  Intra-Conus  models  discussed  in  this  report 
establish  movement  of  forces  to  the  POEs  corresponding  to  these 
chsnnels.  Hence  the  single  piece  of  Information  that  the  Intra-CONUS 
models,  described  in  this  report,  definitely  need  is  allocation  of 
forces  to  POEs. 

Depending  on  the  type  of  Intra-Conus  model  used,  various  other 
pieces  of  information  generated  by  MODES  could  be  used  to  accelerate 
the  solution  process.  In  a  mathematical  programming  based 
decomposition  approach  similar  to  those  discussed  in  earlier  sections 
of  this  report,  certain  additional  information  would  help  in  the 
description  of  t.he  objective  function  or  as  an  advanced  start 
procedure. 

Dual  variables  for  the  time  expanded  channels  would  provide  costs 
associated  with  co-ordinating  the  Intra-Conus  solutions  with  the 
MODES  solutions.  The  duals  for  the  channels  used  would  be  zero  while 
those  for  channels  not  open  would  be  high.  This  cost  Information 
could  be  used  along  with  the  transport  cost  Information  as  the 
objective  of  the  decomposition  models.  Appyling  the  dual  variable 
information  to  the  Intra-CONUS  objective  would  have  the  effect  of 


••king  it  beneficial  for  the  Xntra-CONUS  aodala  to  gravitata  tovarda 
a  solution  similar  to  MODES.  In  fact,  ignoring  any  dlffsrsncss  in 
costs  of  rail  and  sotor  as  vail  as  outload  limitations,  applying 
thasa  dual  valuas  in  tha  Intra-CONliS  objective  of  aithar 
dacomposltion  modal  would  rasult  in  tha  MODES  POE  allocation  baing  an 
optimal  ona. 

Tha  affactlva  tiaa  windows  at  tha  POEs  for  aach  forca,  basad  on 
tha  tima  windows  supplied  by  tha  Supporting  Commander,  could  ba  u sad 
in  tha  objective  function  of  tha  decompositions  models  by  using  costs 
of  zero  for  arriving  within  thasa  windows  and  increasing  costs  for 
arriving  before  or  after  thasa  window  areas. 

An  initial  solution  regarding  tima  of  arrival  of  a  forca  at  its 
destined  POE  could  ba  provided  by  MODES.  This  time-of -arrival 
information  could  ba  used  as  a  starting  solution,  to  ba  improved  upon 
in  subsequent  iterations  of  tha  dacomposltion  models. 

Useful  information  can  ba  provided  by  MODES  to  other  types  of 
Intra-Conus  models  (a.g. ,  deterministic  simulations)  as  follows. 

Tha  time  a  forca  arrives  at  its  POE  in  tha  MODES  solution  could 
ba  used  as  a  target  solution  to  ba  aimed  for  by  tha  modal. 

MODES  might  provide  information  regarding  alternative  POE  choices 
for  a  forca.  This  could  ba  accomplished  by  providing  alternative 
optimal  forca  POE  allocation  scenarios,  or  by  providing  Information 
regarding  those  forca  channel  allocations  that  are  "close”  to  aach 
other. 

Other  types  of  information  would  ba  POE  "equivalence”  relations. 
This  "equivalence"  could  be  established  in  a  varitey  of  ways.  One 
way  would  be  based  on  the  forces  that  MODES  allocates  to  the  POE  over 
time.  Another  way  might  be  by  averaging  the  dual  variables  for  the 
time  expanded  channels  associated  with  the  POE.  This  provides 


additional  flaxibility  to  Intra-Conua  aodaia  by  allowing  (Halted) 


choice  of  a  POE  aaaignaent  which  reault  in  the  leant  tranaport  coat. 
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